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Abstract 

Parton recombination has been found to be an extremely useful model to under- 
stand hadron production at the Relativistic Heavy Ion Collider. It is particularly 
important to explore its connections with hard processes. This article reviews some 
of the aspects of the quark recombination model and places particular emphasis on 
hadron correlations. 
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1 Introduction 

External probes are valuable tools for testing the properties of plasmas. In 
relativistic heavy ion physics hard QCD processes and electromagnetic pro- 
cesses take this role and are are supposed to unravel the secrets of the quark 
gluon plasma (QGP). Indeed the first round of experiments at the Relativis- 
tic Heavy Ion Collider (RHIC) found impressive evidence for jet quenching, a 
dramatic energy loss of the leading jet particle [Tp] . Our interpretation of such 
results relies on the assumptions we make about the applicability of perturba- 
tive quantum Chromodynamics (pQCD). The view that momentum transfers 
of 1 to 2 GeV/c are sufficient is supported by the success of next-to- leading 
order pQCD calculations for p + p collisions at RHIC energies. However the 
answer does depend heavily on the process under consideration. 
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Indeed we have now firm experimental evidence that in Au+Au coUisions at 
RHIC most hadrons are not produced from the usual jet fragmentation mech- 
anism up to rather high transverse momenta Prp ^ 6 GeV/c. This makes 
perturbative hadron production an even rarer process than most of us ex- 
pected. We want to call hadron production perturbative if it originates from 
fragmentation of a jet. Of course jets interact heavily with the surrounding 
medium. In the past it was often assumed that the quenching leads to a reshuf- 
fiing of energy within the jet cone or a transfer of energy from the jet into the 
medium without completely destroying the jet. As it turns out, the number 
of hadrons from surviving jets at intermediate Pt is drowned out by the num- 
ber of hadrons hadronizing from the bulk fireball via recombination of quarks 
[3lill5][6]. 

The hints from experimental data can be divided into three categories. First, 
one-particle inclusive hadron spectra show that the chemical composition of 
the hadrons produced between a Pt of 2 and 6 GeV/c is not compatible with 
jet fragmentation. Most notable is the enhancement of baryons as seen in the 
missing suppression of the nuclear modification factors Raa for baryons and 
in the various anomalously large baryon/meson ratios [7.8. 9|. These results 
already hint toward a chemically equilibrated source of these hadrons. 

Second, the elliptic fiow coefficient V2 follows to very good accuracy a scaling 
law with respect to the number of valence quarks piTUj . This is very direct ev- 
idence that the observable V2 is shaped in a deconfined quark phase. A third 
emerging group of data involves dihadron correlations or associated hadron 
yields. This data is much more difficult to describe (in any model) than the 
single hadron spectra. Although jet-like correlations are present at interme- 
diate transverse momenta, the shape of these correlations and the chemical 
composition can differ strongly from the vacuum [TTP^ . 

The meson has long been discussed as a good test for the validity of the 
recombination model [13j. mesons are as heavy as protons but data from 
RHIC now impressively conffrm that the elliptic ffow and nuclear modiffcation 
factor are very similar to those for kaons We can now be very sure that the 
meson vs baryon signature at intermediate Pt is very robust and can neither 
be explained by an extrapolation of hydrodynamics nor by perturbative jet 
production and fragmentation. 



2 Recombination 

In central heavy ion collisions a hot and dense ffreball of deconffned quarks 
and gluons is created. In the recombination model one postulates the existence 
of thermalized parton degrees of freedom at the phase transition temperature 
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Tc which recombine or coalesce into hadrons. It has been found to be sufficient 
to consider the lowest Fock state in each hadron, the valence quarks, which 
are given constituent masses around 300 MeV. 

The spectrum of hadrons can be calculated starting from a convolution of 
Wigner functions For a meson with valence (anti)quarks a and h we have 

~dFP~^\^J (27r)3 J (27r)3 

xiy„, (^R - ^, I - q; R + ^, ^ + $M(r, q). (1) 

Here Wab is the 2-particle Wigner function for partons a, h and is the 
Wigner function of the meson. The sum runs over all possible parton quantum 
numbers. For simplicity the parton Wigner function is usually approximated 
by a product of single particle phase space distributions Wab = WaWb- Several 
slightly different implementations of this formalism have been discussed in the 
literature j^f5|6] . See [I5] for earlier reviews. 

A very good description of hadron spectra and hadron ratios measured at 
RHIC can be achieved by combining hadron production from recombination 
for intermediate transverse momentum with a perturbative calculation using 
fragmentation and energy loss in the medium [3|H]. Fig. [T] shows the Pt spec- 
trum of 7r°, p, Kq and A + A in central Au+Au collisions obtained in [4j. The 
agreement with data is very good for Pt > 2 GeV/c. We note that the hadron 
spectra exhibit an exponential shape up to about 4 GeV/c for mesons and up 
to about 6 GeV/c for baryons, where recombination of thermal quarks domi- 
nates. Above, the spectra follow a power-law and production is dominated by 
fragmentation. 

Let us now assume the parton phase exhibits elliptic flow V2{pt)- Recombina- 
tion makes a prediction for elliptic flow of any hadron species after recombi- 
nation which is precisely the scaling law mentioned before [TBIH] 

V2{Pt) = nvl{PT/n) (2) 

involving the number n of valence quarks of the hadron. Fig. [2] shows the 
measured elliptic flow f 2 for several hadron species in a plot with scaled axes 
f 2/n vs Py/ n. All data points (with exception of the pions) fall on one universal 
curve. The effect that pions are shifted to lower Pt has been understood [T7] . 
The quark scaling law is an impressive confirmation of the recombination 
model. It shows that the relevant degrees of freedom at early times in the 
collision are partons and it proves that these partons behave collectively. 



3 




Fig. 1. Spectra of tt^, p, Kq and A + A as a function of Pt at midrapidity in 
central Au+Au collisions at ^/s = 200 GeV |4j. Dashed lines are hadrons from 
recombination of the thermal phase, dotted line is pQCD with energy loss, solid 
line is the sum of both contributions. Data are from PHENIX (vr*^, p) [1.7J and 
STAR {K^, A + A) [8]. 
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Fig. 2. Elliptic flow V2 for vr^, p, Kq and A as a function of Pt scaled by the 
number of valence quarks n vs Px/n. The data follows a universal curve, impressively 
confirming the quark scaling law predicted be recombination. Deviations for the 
pions are discussed in the text. Data are taken from PHENIX (vr"^, p) [W\ and 
STAR {K^, A) [9]. 
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Fig. 3. Scaling violation variable A = (B — M)/{B + M) for three different proba- 
bilities 6*2 to have an additional parton beyond the valence quark structure. Wave 
functions of finite width have been used which lead to a scaling violation even in 
the case C2 = of pure valence quark recombination [ISj. 



It is an interesting question to ask to which accuracy one expects the scaling 
law for elliptic flow to hold. In particular, are the scaling factors of 2 and 3 
indeed excluding any higher Fock states in the hadrons? A recent study found 
that higher Fock states in an expansion \p) = aoluud) + ai\uudg) + . . . could 
actually be accommodated [18]. This is in fact trivial for the single-hadron 
yields from a thermal parton spectrum. They do not change if additional 
partons are allowed to coalesce. The situation changes for elliptic flow. Higher 
Fock states with n partons come with their own scaling factor n which seems to 
destroy the scaling with the number of valence quarks. However a numerical 
evaluation shows that the correction is surprisingly small. Fig. [3] shows the 
expected violation of the scaling law using the new asymmetry variable A = 
[B — M) /{B + M) where B and M are the scaled elliptic flow of a meson 
and a baryon respectively [H]. Generally the violations are smaller than 5%. 
New data from STAR analyzes scaling violations in the data and finds them 
to have the predicted sign and order of magnitude [19]. After this study it is 
clear that there is some room to accommodate gluons or sea quarks during the 
recombination process. Further investigations in this direction are necessary. 



3 Hadron Correlations 

One crucial simplification implemented so far is a factorization of any n-parton 
Wigner function into a product of independent single parton distributions 
W^i,...,n = IYi=iWi. By definition this factorization does not permit any corre- 
lations between partons. Consequently, no hadron correlations can emerge via 



recombination. Originally this factorization was chosen for simplicity and it 
was justified because single inclusive hadron spectra could be described very 
well. It has been shown in |5D] that modifications of this factorization includ- 
ing correlations between partons lead to correlations between hadrons after 
recombination. The quality of the single hadron spectra does not suffer. The 
source of jet-like correlations in the QGP is the strong coupling of jets to the 
medium. The energy loss is estimated to be up to 14 GeV/fm for a 10 GeV 
parton [21]. Quenched jets lead to a considerable local heating and create a 
hot spot inside the fireball. Moreover, the directional information of the jet is 
preserved. Partons of such a hot spot exhibit jet-like correlations. 

The next oder extension of the correlation-free factorization in [2Dj assumes 
that correlations still are a small effect and that one can restrict them to 2- 
particle correlations Cij (see also [22]). Then a 4-parton Wigner function can 
be written 

W^i234 ~ ^11^2^3^4(1 + ^ij)- (3) 

i<j 

The correlation functions Cij between parton i and parton j are non- vanishing 
only in a subvolume of the fireball. A Gaussian ansatz Cij ~ cqC"'-'^'"''^^^ /(2</'o) 
is chosen to describe correlations in azimuthal angle. The 2-meson yield is given 
by a convolution of the partonic Wigner function with the Wigner func- 

tions ^A, $B of the mesons with an additional integration over the hadroniza- 
tion hypersurface [20]. It is assumed that the correlation strength cq ^ 1 
which permits omitting quadratic terms like Cq or CqV2- 

We can now study the associated yield Y^b for a trigger hadron A in a given 
kinematic window as a function of the relative azimuthal angle A0 between 
the two. One finds 

27cNaYab{A<^) = gcoe-(^*)'/(20o)^jv^jv^. (4) 

The Ni are single particle yields in the kinematic window of the trigger meson 
or associated meson and cq = cqVc/ (tAt) where tAt is the hadronization 
volume. The factor Q = 4 (for two mesons) indicates an enhancement of the 
correlations in the hadron phase compared to the parton phase. The effect is 
similar to the amplification of elliptic flow by the number n of valence quarks. 
Q counts the number of possible correlated pairs between the (anti) quarks 
of meson A and the hb (anti)quarks of meson B. In the weak correlation limit 
only single correlations are counted. Apparently one has Q = uaUb, thus 
Q = 6 for a meson-baryon pair and Q = 9 for a baryon-baryon pair. 

Fig. m shows the associated yield of hadrons integrated in azimuthal angle 
around the near side (A0 = 0) for the case that the trigger is a baryon (pro- 
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Fig. 4. The associated yield Y^'^^ for baryon triggers (right panel) and meson 
triggers (left panel) as a function of N^art- Squares: fragmentation only; dia- 
monds: fragmentation and recombination with cq = 0; circles: the same with 
Co = 0.08 X lOO/TVpart; data: PHENIX [12]. 

ton or antiproton) and a meson (pion or kaon) for different centralities. The 
kinematic window is 2.5 GeV/c < Pta <4.0 GeV/c for trigger particles and 
1.7 GeV/c < Ptb <2.5 GeV/c for associated particles, and ||/^|, \yB\ < 0.35. 
Fig. m shows the associated yield with only fragmentation, and fragmentation 
and recombination both taken into account together with PHENIX data [12j. 
A good description of the data can be reached assuming a constant correlation 
volume. The parameters used for the fireball are the same that lead to a good 
description of single hadron spectra and elliptic flow [^. 
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